Quantum control of proximal spins using nanoscale magnetic resonance imaging 
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Quantum control of individual spins in condensed matter systems is an emerging field with wide- 
ranging applications in spintronics [1], quantum computation [2] , and sensitive magnetometry [3] . 
Recent experiments have demonstrated the ability to address and manipulate single electron spins 
through either optical [4, 5] or electrical techniques [6-8] . However, it is a challenge to extend 
individual spin control to nanoscale multi-electron systems, as individual spins are often irresolvable 
with existing methods. Here we demonstrate that coherent individual spin control can be achieved 
with few-nm resolution for proximal electron spins by performing single-spin magnetic resonance 
imaging (MRI), which is realized via a scanning magnetic field gradient that is both strong enough 
to achieve nanometric spatial resolution and sufficiently stable for coherent spin manipulations. We 
apply this scanning field-gradient MRI technique to electronic spins in nitrogen- vacancy (NV) centers 
in diamond and achieve nanometric resolution in imaging, characterization, and manipulation of 
individual spins. For NV centers, our results in individual spin control demonstrate an improvement 
of nearly two orders of magnitude in spatial resolution compared to conventional optical diffraction- 
limited techniques. This scanning-field-gradient microscope enables a wide range of applications 
including materials characterization, spin entanglement, and nanoscale magnetometry. 



Magnetic field gradients allow for spatially distinguish- 
ing spins in ensembles, as fields locally modify the spins' 
resonance frequencies. Spatially separated spins can 
therefore be addressed selectively, allowing for magnetic 
resonance imaging (MRI), which has revolutionized the 
medical and biological sciences by yielding micron-scale 
imaging of nuclear spins [9, 10]. Performing MRI on sin- 
gle spins with high spatial resolution is attractive both 
for determining structure on the molecular scale and for 
achieving individual spin quantum control in ensemble 
systems. With conventional MRI techniques, however, 
it is difficult to improve the spatial resolution to the 
nanoscale due to insufficient readout-sensitivity and inad- 
equate magnetic field gradients [11]. Recently, magnetic 
field gradients introduced via scanning probe techniques 
have enabled single spin detection with few-nm resolu- 
tion [12, 13]; however, control and characterization of in- 
dividual spins in nanoscale clusters has not been demon- 
strated thus far. 

Here we perform scanning field-gradient MRI on prox- 
imal electron spins in nanoscale ensembles and demon- 
strate a spatial resolution < 10 nm under ambient condi- 
tions. We show that scanning field- gradient microscopy 
not only allows for imaging but further provides quan- 
tum spin control for characterization and manipulation 
of individual spins on the nanoscale. By pushing the 
spatial resolution to few-nm length scales, our results il- 
lustrate that quantum control of individual spins can be 
maintained in dense ensembles of spins, where the mu- 
tual coupling between adjacent spins can become very 
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strong [14]. Thus, scanning- field-gradient MRI will help 
facilitate the creation of entangled spin-states with ap- 
plications for quantum information processing and sensi- 
tive, nanoscale magnetometry. While our approach is ap- 
plicable to any spin system where spins can be initialized 
and read out, we focus here on the electronic spins asso- 
ciated with NV centers in diamond, where spin initializa- 
tion and readout can be performed optically [15]. Addi- 
tionally, NV spins are attractive for performing quantum 
information processing [14, 16] and sensitive magnetome- 
try [13, 17-19]. Individual NV spin control in a nanoscale 
ensemble is a key advance towards the implementation of 
these applications. 

Our scanning- field-gradient MRI system (Figure la) is 
comprised of an atomic force microscope (AFM) outfit- 
ted with a magnetic tip and integrated into an optical 
confocal microscope, all operating under ambient condi- 
tions. Small ensembles of shallowly implanted NV cen- 
ters, ^10 — 50 nm below the surface of a diamond sam- 
ple, are placed in the confocal spot (volume < 1 cubic 
micron) where an excitation laser at 532 nm is used 
to initialize and read out the NV centers' spin states. 
Because of the NV centers' spin-dependent fluorescence 
(~ 630 — 800 nm), optically detected electron spin reso- 
nance (ESR) can be observed by sweeping the frequency 
of a driving radio-frequency (RF) field through the spin 
resonance and measuring the corresponding variation in 
fluorescence (Figure lb) [15], in this case on a single NV 
center. As has been demonstrated by Balasubramanian 
et al. [13], a magnetic tip in proximity to an NV center 
shifts the energy of the NV spins, particularly for fields 
along the NV axis. By selectively detuning the applied 
RF field and scanning the magnetic tip over a single NV 
center, a magnetic field map corresponding to the applied 
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FIG. 1: Two-dimensional imaging of a single NV center using 
the scanning-field-gradient MRI microscope, (a) Schematic of 
the experimental approach. A magnetic tip is scanned by an 
AFM over a diamond sample containing multiple shallow NV 
centers separated by distances smaller than the optical excita- 
tion wavelength. Each NV center experiences a different mag- 
netic field from the tip. Thus, spin transitions for individual 
NV centers can be selectively driven by tuning the frequency 
of an externally applied radio-frequency (RF) field. A confo- 
cal microscope provides NV spin-state preparation via optical 
pumping, and detection via spin-state dependent fluorescence, 
(b) Optical fluorescence measurement of the electron spin res- 
onance (ESR) of a single NV center in the absence of the mag- 
netic tip. The NV-spin-dependent fluorescence rate leads to 
a drop in emitted (and detected) photons when the external 
RF source is swept onto the NV spin resonance (here at 2.47 
GHz instead of the zero-field value of 2.87 GHz because of an 
applied static magnetic field), (c) Two-dimensional magnetic 
resonance image of a single NV center, created by scanning 
the magnetic tip across the surface and fixing the RF fre- 
quency off resonance from the NV ESR transition (see arrow 
in (b)). A reduction of fluorescence is observed for positions 
of the magnetic tip relative to the NV that put the NV spin 
on resonance with the RF field, creating a dark "resonance 
ring". The plotted fluorescence magnitude is normalized to 
the NV fluorescence when no driving RF field is applied near 
the ESR transition. 



detuning can be acquired (Figure Ic). Using the AFM 
to regulate the magnetic tip's height allows for precise 
nanoscale control of the tip's location in three dimen- 
sions: the height of the tip is maintained via a conven- 
tional feedback loop on the force from the sample, while 
several-nm precision in lateral positioning is achieved by 
using the sample's topography as a marker. 

We first demonstrate how this technique can be used 
for imaging proximal NV spins. The magnetic response 



shown in Figure Ic provides a direct means for determin- 
ing the relative location of proximal NV centers. For 
NVs with an identical orientation, the indistinguishable 
nature of their spins leads to an identical but spatially 
shifted magnetic field map for each NV spin. Thus, mag- 
netic resonance images of single-spins (such as in Fig- 
ure Ic) can serve as the point spread function for multi- 
spin imaging. To demonstrate the performance of this 
technique, we executed scanning- field-gradient imaging 
on three closely spaced NV centers (Figure 2). As the 
magnetic tip is scanned laterally across the sample, a 
magnetic field contour is observed for each center. Here, 
we selected an ensemble of NVs where all three have a 
common quantization axis orientation, so that every spin 
responds in the same way to the presence of the magnetic 
tip. Thus, the relative distances between NV spins can 
be obtained by quantifying the spatial shift between the 
observed resonance rings through a deconvolution proce- 
dure. The resulting image (Figure 2a) indicates the rela- 
tive positions of the NV centers, which we find are spaced 
by 50 nm and 75 nm with respect to NV H. A spatial 
resolution of 9 nm can be extracted from the width of 
the resonance ring along the vector connecting NV cen- 
ters n and HI; a precision of 0.2 nm can be determined 
from a two-dimensional Gaussian fit to the deconvolved 
peak. This estimate of precision only accounts for vari- 
ance induced by random noise and does not account for 
systematic deviations caused by effects such as scanner 
non-linearity or deviations from the peak shape to the 
Gaussian fit. 

Our single-spin MRI technique can be extended to 
three dimensions by performing magnetic tomography, 
where the magnetic tip is retracted from the surface in 
few-nm steps and scanned laterally (Figure 2b) . The size 
of the resonance rings evolves quickly as a function of 
z-distance as the magnetic field from the tip becomes 
too weak to bring the spin-transitions into magnetic res- 
onance with the detuned RF driving field. The height 
differences between the NVs can be seen here as NV I 
vanishes roughly 15 nm before the other two, which in- 
dicates that NV I lies roughly 15 nm further below the 
surface than both NV II and III. A vertical scan across 
NVs II and II shows the z-resolution of this measure- 
ment to be 10 nm, which is extracted from the width of 
the resonance line in the z-direction. Since the lateral 
separation of these two NVs is larger than our spatial 
resolution, we can determine the relative height of the 
two NVs with high precision, which we find to be 3 nm. 

The NV centers in our demonstration experiments 
were created through implantation of nitrogen ions form- 
ing a layer ^10 nm below the diamond surface. Model- 
ing of this implantation procedure [20] predicts a spatial 
variation of ±3 nm; whereas we observe a larger varia- 
tion in the distribution of NV depths (> 10 nm) using 
our scanning- field-gradient MRI technique. This discrep- 
ancy may arise from ion channeling or surface effects. 
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FIG. 2: Three-dimensional magnetic resonance imaging of proximal NV spins, (a) Scanning the magnetic tip over a cluster of 
three NV centers with the same crystallographic orientation yields multiple dark-resonance rings in the observed fluorescence, 
one for each NV center (upper image). The relative locations of the NV centers are extracted from the magnetic resonance image 
through a deconvolution process, yielding adjacent NV-NV distances of 50 nm and 70 nm (lower image). ESR spectral linecuts 
(e.g., along x' as shown here) give a spatial resolution of roughly 9 nm with 0.2 nm precision taken from a two-dimensional 
Gaussian fit to the deconvolved peak, (b) The relative depths of the NV centers below the diamond surface (i.e., relative NV 
positions along the z-axis) are determined by taking magnetic resonance images for different magnetic tip heights above the 
sample surface. Comparing the evolution of the three NV resonance rings as a function of tip-to-sample distance, we determine 
that NV I lies roughly 15 nm below NV II and NV III (stack of images on left). A vertical cut further resolves NV II and NV 
III, showing that they are 3 nm apart in depth (right image). 



Nanoscale precision in measuring the distance between 
proximal NV spins also allows determination of the mu- 
tual dipole coupling between adjacent spins, a key com- 
ponent for creating entangled spin-states. 

Scanning-field-gradient MRI provides not only a 
method for nanoscale mapping of the spatial locations of 
proximal spins but also allows for individual spin tran- 
sitions to be resolved in frequency space with high pre- 
cision. In the absence of magnetic field gradients, iden- 
tical spins sharing a quantization axis are indistinguish- 
able making selective control over proximal spins impossi- 
ble. Performing scanning- field- gradient MRI on proximal 
spins differentiates their transition frequencies, allowing 
for coherent manipulation and characterization of indi- 
vidual spins in the ensemble. This can be done simulta- 
neously on several NV spins while preserving independent 
control of each spin. Such selective control over proximal 
spins can be maintained as long as the magnetic field 
gradient separates the transition frequencies of neighbor- 
ing spins by more than their resonance linewidth, which 
is the same condition that determines the spatial reso- 
lution in MRI. Therefore, spins separated by more than 
the achieved MRI spatial resolution (^ 9 nm for the NV 
spin experiments reported here) can be addressed inde- 



pendently using our technique. 

To demonstrate such selective nanoscale characteriza- 
tion and control of proximal spins we examined a pair 
of proximal NV centers (NV IV and NV V, separated by 
135 nm, which share the same NV axis orientation. In the 
absence of the magnetic tip, we performed a continuous- 
wave ESR measurement on both NV spins simultane- 
ously (Figure 3a, left panel), finding no difference in their 
spectra, as expected. After tuning the RF frequency to 
the pair of NV spins, we drove Rabi oscillations and ob- 
served the free induction decay of the two NV spins via 
a Ramsey sequence (Fig 3a, right panels). The measured 
Ramsey fringes show a pronounced beating pattern due 
to the hyperfine structure of the NV spin transitions and 
are damped due to inhomogeneous dephasing of the two 
NV centers. However, information distinguishing the in- 
dividual coherence-properties of the proximal NV centers 
cannot be extracted from these measurements. In con- 
trast, when the magnetic tip is brought in close proxim- 
ity to the NV spins, it splits their resonance frequencies 
due to the differing magnetic fields applied to each NV, 
thereby allowing each NV spin to be addressed and char- 
acterized individually and coherently (Figure 3b) without 
modifying the spin of the neighboring NV. 
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FIG. 3: Individual quantum control and characterization of 
proximal NV centers, (a) With the magnetic tip pulled far 
away from the sample, two NV centers (IV and V) separated 
by 135 nm and sharing the same spin-quantization axis cannot 
be resolved by their ESR spectra (ms = ^ ms = 1 tran- 
sition of the NV triplet groundstate, observed by spin-state 
dependent fluorescence), as they experience the same static 
magnetic field (left panel). By driving both NVs at once, 
Rabi oscillations and the NV spins' collective free-induction 
decay can be observed using the indicated Rabi and Ram- 
sey RF pulse sequences, respectively; but the evolution of 
each individual NV-spin remains undetermined, (b) With 
the magnetic tip in close approach to the sample, the two NV 
spin resonances are spectrally distinguishable by their differ- 
ing Zeeman shifts (~ 30 MHz splitting), allowing each NV 
to be addressed independently, so that Rabi oscillations and 
Ramsey free-induction decay of each NV spin can be indi- 
vidually measured (NV IV in red and NV V in green). For 
these measurements, the probability of the spin being in the 
ms — state is plotted. We observe that NV IV has a shorter 
inhomogeneous coherence time (T2 ~ 180 ns) than NV V 
(T2* ~ 420 ns). The Ramsey free induction decay measure- 
ments are taken with the RF driving field detuned by 5 MHz 
from the target NV center's nominal ESR frequency in the 
presence of the magnetic tip; the observed oscillations are due 
to beating of this detuning with the 15N hyperfine splitting 
(3.1 MHz). The measured Ramsey data are fit to the sum of 
two exponentially damped sinusoids, whose phases are fixed 
by the relative strength between the net detuning (the RF 
field detuning plus or minus half the hyperfine splitting) and 
the Rabi frequency (5.5 MHz), governed by the strength of 
the applied RF field. 



Such coherent individual spin control was realized by 
tuning the RF frequency to the ESR resonance of the tar- 
get NV center. Pulsing the RF with a variable duration 
induced coherent Rabi oscillations of either NV IV or NV 
V, dependent on the tuning of the RF frequency. Simi- 
larly, characterization of individual spin coherence prop- 
erties was achieved by measuring the spin's free-induction 
decay via a Ramsey sequence using the appropriate RF 
frequency for the target NV spin. We observed that NV 
IV has a faster free-induction decay rate than NV V, indi- 
cating that the two spins have different inhomogeneous 
dephasing times, T2*. The measured free-induction de- 
cay rate for NV V with and without the magnetic tip 
are comparable, indicating that additional decoherence 
induced by the tip is small compared to that spin's am- 
bient dephasing 1/T|. In order to prevent tip-induced 
spin-decoherence, we employed a three-dimensional spa- 
tial feedback scheme, which ensured that variations in the 
applied tip-field were smaller than the intrinsic NV inho- 
mogeneous dephasing rates. Thus, the applied magnetic 
field gradient can be used to characterize individual spin 
coherence properties, a direct consequence of achieving 
single-spin control. 

Under the influence of the tip's magnetic field gra- 
dient, manipulating one NV spin does not perturb the 
state of a neighboring NV spin. To experimentally verify 
this selective, independent spin control, we simultane- 
ously drove Rabi oscillations on NV IV while measuring 
the ESR spectrum of NV V (Figure 4a) by applying two 
separate RF fields (RF 1 and RF 2), each used to drive 
one NV spin-resonance. To illustrate the independence of 
the two measurements, we performed a two-dimensional 
sweep over the RF 1 pulse duration {trfi) and the fre- 
quency of RF 2 {(jjRF2) (Figure 4b). We then subtracted 
out the mean of each row or column to show that the two 
measurements are independent. Rows or columns were 
summed to reconstruct the resulting Rabi oscillations of 
NV IV or the ESR curve of NV V, respectively. The only 
deviations from simultaneous, independent NV spin mea- 
surements are observed when (jJrf2 is swept close to the 
resonance of NV IV, yielding a slight damping in the Rabi 
oscillations. These deviations do not persist through the 
ESR transition of NV V, which shows that both NVs can 
be manipulated independently when they are driven on 
resonance. 

Scanning-field-gradient MRI of spins yields precise de- 
termination of their relative locations with nanomet- 
ric spatial resolution, which will be crucial for optimiz- 
ing the performance of spin-based magnetometers and 
the functionality of spin-based quantum bit ensembles. 
Once suitable spin ensembles are identified and spatially 
mapped, individual spin control both allows for the deter- 
mination of individual spin properties and - when com- 
bined with magnetic dipole coupling between adjacent 
spins - provides a method for achieving complete con- 
trol of the quantum state of spin ensembles. The spa- 
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FIG. 4: Selective, independent RF control of proximal NV 
spins in the presence of the tip magnetic field gradient, (a) 
Fluorescence measurements are performed for the same two 
proximal NV centers as in Figure 3, while undergoing simul- 
taneous, near-resonant driving by fields RF 1 and RF 2; 
and with the tip magnetic field gradient inducing a large 
30 MHz) Zeeman frequency shift between the ESR fre- 
quencies of the two NV spins. The frequency of RF 1 is set 
on resonance with NV IV and is pulsed with varying dura- 
tion (trfi) to induce Rabi oscillations. Simultaneously, RF 
2 is continuously applied while its frequency {ujrf2) is swept 
through the spin resonance of NV V to measure its ESR spec- 
trum, (b) To determine if these simultaneous NV measure- 
ments are independent, a two-dimensional data set containing 
all relevant values of both trfi and ujrf2 is acquired (base of 
cube). From this data set, the individual behavior of each NV 
can be extracted by subtracting the mean measured fluores- 
cence for each row (i.e., fixed value of ujrf2) from the data for 
that row; and similarly for each column (fixed value of trfi). 
The resulting extracted data sets (shown in the sides of the 
cube) are found to be independent of row or column number, 
showing that varying the Rabi pulse duration trfi for NV IV 
does not influence the results of sweeping ujrf2 through the 
spin resonance of NV V, and vice versa. The only observed 
deviations from this independent NV spin control are for small 
values of ojrf2 where the frequency of RF 2 approaches the 
spin resonance of NV IV, thereby causing modest damping of 
this NV's Rabi oscillations. These deviations disappear near 
the resonance of NV V, showing that both NVs can be driven 
independently. Summing up the measurements for all rows 
and columns yields the resulting averaged Rabi oscillation 
measurement for NV IV (top of left wall) and ESR spectrum 
from NV V (top of right wall), respectively. 



tial resolution of our scanning technique is extendable to 
the atomic scale by using stronger magnetic field gradi- 
ents [21] and narrower ESR linewidths [22]. For NV spins, 
the current experimental barrier to improving this spatial 
resolution is overcoming the reduction in ESR contrast 
due to a strong off-axis magnetic field created by the 
tip [23]; however, this can be prevented by either apply- 
ing a strong bias-field along the NV axis or by tailoring 
the domain structure of the magnetic tip to produce high 
gradients but only moderate total fields. 

The control and manipulation of individual spins 
using magnetic field gradients is independent of the 
method used for spin readout. For optically address- 
able spins, such as NV spins, integrating far-field, sub- 
diffraction schemes - such as stimulated emission de- 
pletion (STED) [24] and reversible saturable optical lin- 
ear fluorescence (spin-RESOLFT) [25] - with a scanning 
magnetic field gradient would allow for both robust indi- 
vidual spin control and readout with nanometric resolu- 
tions. Additionally, selective optical control of such sys- 
tems is possible via the incorporation of an electric field 
gradient to the scanning tip [26] , which would allow both 
spin and electronic degrees of freedom to be both ad- 
dressed individually. Alternatively, using demonstrated 
single-shot electrical readout of individual spins [7, 27] 
would allow for MRI to be performed rapidly and effi- 
ciently, as acquisition times would not be limited by the 
readout integration time. Individual control of spins via 
the present technique is also extendable to nuclear spins, 
provided methods for reading out nuclear spins reach 
single-spin sensitivity [21]. Nanoscale or atomic spatial 
resolution of nuclear spins is feasible as their long spin co- 
herence times help to compensate for nuclear spins' small 
dipole moment. For any spin system, providing individ- 
ual control of spins in dense ensembles, where mutual 
coupling is strong, allows for the creation of arbitrary en- 
tangled states. Such states have intriguing potential ap- 
plications ranging from sensitive nanoscale magnetome- 
ters to scalable quantum information processors [28]. 
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Methods 
NV center samples: 

NV centers were created through implantation of 15N 
ions [29] into ultrapure diamond (Element Six, electronic 
grade diamond, < 5 ppb nitrogen). The implantation 
was done at 6 keV to give a nominal nitrogen depth of 
10 nm [20]. To form NV centers, the sample was annealed 
in vacuum at 750 °C, where existing vacancies are mobile 
and can pair with the implanted nitrogen atoms. The re- 
sulting density of NV centers corresponds to one center 
every 50-100 nanometers, forming a layer roughly 10 nm 
from the surface. To isolate small NV clusters (i.e., a 
few proximal NV centers) , we selectively etched [30] away 
the majority of the shallow diamond surface layer, leav- 
ing proximal NV-containing nanostructures. This was 
done using electron beam lithography to define an etch 
mask from a flowable oxide (Dow Corning, XR-1541) [31]. 
A reactive-ion etch then removed any exposed diamond 
surfaces, resulting in shallow diamond nanostructures 
(100 — 800 nm across) containing ensembles of proximal 
NV spins. 

Magnetic tips: 

Magnetic tips were created by evaporating a magnetic 
layer onto quartz tips of roughly 80 nm in diameter, 
which were fabricated using a commercial laser-pulling 
system (Sutter Instrument Co., P-2000). Using a thermal 
evaporator, a 25 nm layer of cobalt-iron was deposited on 
the side of the pulled quartz tip. A 5 nm chrome layer 
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was then evaporated, which serves as a capping layer to 
prevent oxidation of the magnetic material. These tips 
result in magnetic field gradients of roughly IG/nm at 
distances of roughly 100 nm. 

Magnetic tip positioning with an AFM: 

To achieve high spatial resolution in NV imaging and 
manipulation, it is necessary to control precisely the rela- 
tive distance between the magnetic tip and the addressed 
NV centers. Such nanoscale control can be challeng- 
ing under ambient conditions (standard temperature and 
pressure); e.g. temperature drifts of a small fraction of 



a degree can induce few nm drifts between the magnetic 
tip and the diamond sample, which would inhibit local 
spin-control. 

To overcome this problem, we used an AFM to posi- 
tion the magnetic tip in three dimensions with precisions 
of a fraction of a nanometer in z and a few nanometers in 
both X and y. Height control was achieved through nor- 
mal AFM operation in which a feedback loop modulates 
the height of the tip to keep the sample-tip interaction 
constant. Lateral positioning was achieved through in- 
termittently locating topographic features on the sample 
to ensure that the tip's relative distance to the sample is 
fixed during any performed measurements. 



